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Abstract 

We investigate the effective interaction between Cooper-pair molecules in the strong-coupling 
BEC regime of a superfluid Fermi gas with a Feshbach resonance. Our work uses a path integral 
formulation and a renormalization group (RG) analysis of fluctuations in a single-channel model. 
We show that a physical cutoff energy oj c originating from the finite molecular binding energy 
is the key to understanding the interaction between molecules in the BEC regime. Our work 
thus clarifies recent results by showing that om = 2of is a bare molecular scattering length while 
a-M = (0.6 ~ 0.75)of is the low energy molecular scattering length renormalized to include high- 
energy scattering up to lo c (here of is the scattering length between Fermi atoms). We also include 
many-body effects at finite temperatures. We find that om is strongly dependent on temperature, 
vanishing at T c , consistent with the earlier Bose gas results of Bijlsma and Stoof. 

PACS numbers: 03.75.Ss, 03.75.-b, 03.75.Kk 
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Recently, the interaction between bound molecules has attracted much attention in the 
strong-coupling BEC regime of a superfluid Fermi gas DDDD In a trapped Fermi gas 
with a Feshbach resonance, one can tune the magnitude of a pairing interaction by varying 
the threshold energy of the Feshbach resonance [a, 0, Using this tunable interaction, the 
BCS-BEC crossover|7J has been observed in 40 K and 6 Li, where the character of superflu- 
idity continuously changes from the weak-coupling BCS type to the BEC of tightly bound 
molecules, as one increases the pairing interaction [l], 8^ In the strong-coupling BEC regime, 
where the molecules already form above the superfluid phase transition temperature T c 0], 
the system can be regarded molecular Bose eas. 

The magnitude of the molecular interaction in the BEC regime was originally studied 
in the superconductivity literature j]} Q], where the result a M = 2a F was obtained using a 
Ginzburg-Landau (GL) type expansion. Here o-m is the molecular scattering length and 
ap is the atomic s-wave two-body scattering length. Later, Pieri and Strinati[3] obtained 
<2m — 0.75of within a two-body t-matrix approximation. More recently, Petrov and co- 
worked obtained a M K 0.6a F by solving a fonr-fennion problem. This last resn.t is 
confirmed by a direct Monte-Carlo simulation for the ground state[Q(. In addition, Grimm 
and co-workers also obtained aM — 0.6ap in superfluid 6 LiQ], using a Thomas- Fermi fit to 
their molecular condensate density profile. 

In this paper, we investigate this interaction between molecules in the strong-coupling 
BEC regime of a uniform Fermi gas, using a single-channel model. We show that the finite 
binding energy of a Cooper-pair "molecule" naturally leads to a physical cutoff energy u> c , 
which is not present in the usual atomic BEC. In a two-molecule case, we show that u c is 
the key to understanding the physics behind the difference between earlier results, namely, 
gsm = 2ap |li Q] and au = (0.6 ~ O.75)aF0, 4\. We also calculate the many-body t-matrix 



at finite temperatures, extending the renormalization group analysis for an atomic BEC 
to the BEC regime of a superfluid Fermi gas. In a many-particle system, we show that om 
depends strongly on temperature, vanishing in the region near T c . 

We consider a two-component Fermi gas described by the single- channel BCS model. 
In the functional integral formalism, the action S in the partition function Z = 
£ ff JZ^tp^e-s has the form | 2 | j S = fgdrfdr^ *t (jL + g _ ^ _ U^\%^' . 

Here, ^(r, r) and ^(r, r) are a Grassmann variable and its conjugate, describing Fermi 
atoms with pseudo-spin o =f, j. fi is the atomic chemical potential. In 40 K and 6 Li Fermi 
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gases, the tunable pairing interaction U is associated with the Feshbach resonance 

m. 

We simply treat U as a tunable parameter in a single channel formulation. This is valid for 
a broad Feshbach resonance. 

To discuss the strong-coupling BEC regime, it is convenient to introduce a Cooper-pair 
Bose field A(r, r), using the usual Stratonovich- Hubbard transformation Q]. After functional 
integrations over ty a and we obtain Z = j T>/S)T>/S.e~ SA , where Sa = Jo dr J dr^j 



Tr 



ln[— G x ] [2J. The fermion single-particle Green's function is given by 




(1) 



where Tj (j = 1, 2, 3) is the Pauli matrices. The mean-field gap equation is obtained from 
the saddle point solution for A M f determined by dS^/dA = 0. Expanding the action Sa 
around this mean-field solution (Amf), we obtain 
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where SA(r, r) = A(r, r) — Amf describes fluctuations in the particle-particle channel. The 
mean-field 2x2 matrix Green's function G° has the standard BCS form, 



G°(p,iu ri 



(e p -fi)r 3 - A MF n 



(3) 



Here, iu m is the fermion Matsubara frequency. E p = 



fi) 2 + A MF is the single-particle 



excitation spectrum, where e p is the kinetic energy of a Fermi atom. 

In the BEC regime, the atomic chemicalpotential \x is related to the atomic scattering 
length ap by the formula p, = — l/2ma|Q, 2|- ^ e no ^ e ^ na ^ ft * s l ar g e an d negative as we 
approach the BEC limit (a F — > +0). In this regime, since the BCS order parameter at T = 
is given by Amf = \J16/3tc Ift^^ep 4 ^, we find |/2| ^> A M f- Thus, we can use approximation 
£ p + in (PJ). In addition, the energy gap E g is given by £J S = J |/i| 2 + A MF ~ |/2|, 



~ 



so that the binding energy -Ebind of a Cooper-pair is given by -Ebind = 2E g = 2|/i|. 

In the BEC regime, we expand (j2J) inpowers of 8A, retaining terms up to n = 4. We 



recall that the Gaussian approximation 




only keeps fluctuation terms to order 
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(a) V ( ^(q,q',K) 



(b) Vj^K) 



FIG. 1: (a) Interaction between molecules VL^ mediated by Fermi atoms. The dotted lines 
describe the molecular Boson propagator, (b) Particle-non-conserving interaction Vjjr ■ 



n = 2. In the superfluid phase, we find 
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Here, Um = 47r(2ap)/M (where M = 2m) is the bare s-wave molecular interaction obtained 



in what amounts to the time-dependent GL theory 



q 2 /2M is the kinetic energy 



of a free molecule, and zz/ n is the pair boson Matsubara frequency, with p = (p, iu m ). 
In the BEC limit of interest, the fluctuations in the order parameter SA q involve the 
molecular condensate. We introduce a renormalized molecular field $ 9 = 5A q /i], where 
i] = [X) p l/4(e p — /i) 2 ] -1 / 2 . Defining = Amf/v as the equilibrium value of the Bose- 
condensate order parameter, at T = 0, one can show that n c = <p 2 M = rip/2 (where hf is the 
number of Fermi atoms) gives the number of Bose-condensed molecules. 

V^g in (jl]) describes an interaction between molecules mediated by the Fermi gas [see 
Fig. 1(a)], given by 



Q-q' K Q q + q' K s 

— + 4-) G n(-P+ — + 4", 



4 



1 2 3 4 5 

qa F 

FIG. 2: Calculated momentum dependence of the boson-boson interaction vf$ (q, 0, 0) (solid line). 
The dashed line is the Lorentzian fit given by Um/0- + (<z/A) 2 ). The vertical dotted line shows the 
position of the cutoff momentum A. 

xG^-p-i^ + fx^-^ + f). (5) 

Here, G^i is the diagonal (l,l)-component of the BCS mean-field Green's function defined 
in (jSJ). V}g in (JJJ) represents a particle-non-conserving interaction, given by [see also Fig. 
1(b)] 

V$\q, iO = "jE G°n(P + j)G° u (-p + j)G° 12 (p + q). (6) 
This interaction involves the off-diagonal (anomalous) Green's function Gi 2 denned in (J3j), 

(3) 

so that Vg ff only exists in the superfluid phase. 

For q = q' = K = 0, one can show that ^(0,0,0) = U M and Vj?(0,0) = y/W c U M , 
where Um = 47r(2ap)/M. Using these results in place of V^\q, q', K) and v}jp(q, K) in 
(0J), we reproduce the previous results obtained in the literature, namely, om = 2ap 
However, v}q and actually depend on momentum, and go to zero with a characteristic 
cutoff momentum A. To calculate this cutoff momentum A, we expand (0) and (EJ) in terms 
of the relative momentum (taking the center-of-mass- momentum K to be zero 12|), keeping 
terms of 0{q 2 ). We obtain V$(q, q', 0) ~ U M F{q)F{q') and V}g\q, 0) ~ ^/n~ c U M F{q), where 
F(q) = 1 — (g/A) 2 . The characteristic momentum is found to be A = 4/\/5aF. For larger 
values of the momentum, we replace F(q) by the Lorentzian [1 + (q/A) 2 }^ 1 . This gives a very 
good fit to the momentum dependent interaction, as shown in Fig. 2. In terms of energy, we 



can write F(q) = [1 + /uj c ] 1 , where the interaction or molecular cutoff energy is defined 
by 

A 2 8. . 

" c= 2M = 5 l/i| - (?) 
We note that u> c is comparable to the binding energy Eund = 2|/z| of a molecule in the BEC 
limit. Thus we conclude that u c is a physical cutoff determined by the dissociation energy 
of a molecule in the BEC phase. 

To summarize, when we replace V}^(q, q', K) and v}jp(q, K) in (@J with IQj (0, 0, 0) and 
V}g (0, 0), respectively, we have to include the physical momentum cutoff A, as discussed 
above. We recall that $ g = SA q /i] describes excitations from the Bose-condensate. The 
molecular Bose field involving the condensate part is given by <p q = $> q + <pM- Using (f> q in 
(JUJ), we obtain 

P q 

1 4ir(2a F ) A , t , + 

(3) 

where /im = Umk-c is the molecular chemical potential. We note that V^ s has been absorbed 
into the last term of (jHJ. The corresponding Hamiltonian has the form 

#BEC = H^q( £ f _ ^A/)0q 

q 



47r(2a F ) 



A 



E 4>% +al <t>%_ a ,^^ q (j)K + (9) 



re two different molec- 
and om = (0.6 ~ 



M T+ 9 ' a" 9 ' t_ * 

This Hamiltonian is also obtained in the normal phase, but now with Hm — 4(/a + 

The effective Hamiltonian in gives a clear physical picture of t 
ular scattering lengths found in earlier literature, namely om = 2ap 

O.75)aF0, Q]- Although Um = 47r(2ap)/M looks like the low-energy expression for s-wave 
interaction, (JSJ) shows that Um is actually a 6are interaction involving a cutoff energy u; c . In 
the case of two molecules in a vacuum, to obtain the low-energy interaction Uy = 4Tra^/2M, 
we have to renormalize Um to remove the need for a high-energy cutoff In the usual 
manner, this is given by 

4tt< _ 4tt 2a F 



M i + ±fe) E - c _J^_ Ml + 2.28 

= 47r(0.61a F ) 

M 1 ' 



Here we note that ap in the denominator has been canceled out by the of involved in u c , 
defined in (|7jl. The resulting renormalized Hamiltonian has the same form as 0, but now 
with 0.61of in place of 2a-p 14]. Eq. (jl(Jj) agrees with the recent four-fermion (two-molecule) 
analysis, which gives au = 0.6ap 4|. Thus, we find that au = 2ap is the bare molecular 
scattering length before renormalization, while cim = (0.6 ~ 0.75)ap is the renormalized 
two-molecule scattering length in which the high-energy processes up to u c have been incor- 
porated. 

The effective action in (jHJ) is also useful to calculate the molecular scattering length om 
in a many-particle system. In this case, clm is affected by many-body effects, as well as by 
temperature, which are not included in the T = two-molecule result in (jl(J|) . To include 
these effects, we apply the RG theory developed for an atomic Bose-condensed gas 1]| to 
deal with®, including the physical cutoff A. In the 1-loop level, the RG equations are 
given by [11], for fi M < 0, 
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7T Z 



dl 



~—U M N{sf 



(11) 



dU, 



M 



dl 



-u 



M 



A 3 

—U 2 
2^ M 



r2N(e 



M 
A 



fJ'M) 



+ Af3N{ef -ti M )[N(e 



2(ef 

M 



Vm) + 1] 



and for fiM > 



Q, 



dn 



M 



dl 



2/iM — 



2tt 2 

M 



u 



M 



[2AT«) + 1] 



+ WN{E%)[N{E%) + \] 



(12) 



(13) 



dU, 



M 



dl 



-u 



M 



A 3 

—II 2 

O 2 M 

2n z 



2N{Ef) + 1 



2Ef 



+ ^N{Ef)\N{Ef) + \\ 



Here, Ef 



M( C M 



(14) 



,£ A [ £ a + 2/iA/), and N(e) is the Bose distribution function. Um, and T 
involve the trivial scaling as Hm(1) = Hm^ 21 , Um(1) = Um^~ 1 , T{1) = Te 21 . We solve the RG 
equations together with the equations for the number of molecules hm = uf/2 = n c + n n . 
Here, n c = lim^ oo /i M (0/^ r M(0 e_3 ' (a*m(0 > 0) is the number of Bose-condensed molecules, 
and the number of non-condensate molecules n n is determined bvjll|. 

A 3 



n. 



2tt 2 Jo 



dlN(ef - /iA/)e~ 3 ' (jjl m < 0) 



(15) 
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FIG. 3: Many-body molecular scattering length clm as a function of temperature, p-p is the Fermi 
momentum for a free Fermi gas. We show results down to T = 0, although the RG equations ()13|) 
and (|14j) are not valid for temperatures such that T/Tp <, 0.1(pF fl F)-The inset shows T c in the 
BEC regime determined from the fixed points of the RG equations (|13|) and (|14|) . The dashed line 
shows T c ideal = 0.218T F . 
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Figure 3 shows the many-body molecular scattering length au at finite temperatures. The 
scattering length a M depends on temperature, and the temperature dependence is stronger 
as one approaches the crossover regime. At T c , om vanishes. This is expected since if 
molecules are to be stable at T c , the scattering between them must vanish at T c . We note 
that a similar T-dependent scattering length was also obtained earlier in the study of Bose- 
condensed gases jll|, |l6j. However, in the present situation, one may have such large values 
of dp near the resonance that the decrease in the effective value of cim may have observable 
consequences at finite temperatures. In atomic Bose gases, the small size of the scattering 
length meant this region was restricted to temperatures very close to T c and hence of less 
interest. 

As shown in the inset in Fig. 3, the many-body effects enhance T c compared with T c ldeal = 
0.218Tp for the transition temperature of a non-interacting Bose gas. In the extreme BEC 
limit, we find AT C = T c — T c ldcal oc Qof^f) oc (n^aM), consistent with the previous work on 
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the BEC in Bose gases Il7l Il8l. Il9l| . Although T c in the inset initially increases with p^a? 
from the BEC side, it will eventually decrease in the crossover region, smoothly going into 
the weak-coupling BCS result [3]. 

Our theory is only valid in so far as the effect of fluctuations involving terms higher order 
than n = 4 in (j2J) is small. Comparing the fourth order term with the sixth order term 
above T c , we find that the fourth order term is dominant as long as (pfOf)" 1 ^ 1, i-e., we 
are not in the unitarity region. 

To conclude, we have studied the effective interaction between bound molecules in the 
BEC regime of a superfluid Fermi gas. We have included the binding energy of a molecule, 
which naturally leads to a physical cutoff energy in the molecular Bose gas. For the two- 
molecule case, we showed that this cutoff gives a simple physical explanation of the difference 



between the two molecular scattering 



engths, the bare value ayi = 2of 
malized value czm = (0.6 ~ 0.75)ap|3|, |4j. We have also shown that the many- body scattering 
length deviates considerably from the two-body scattering length as the temperature ap- 
proaches T c if we are close to the unitarity limit. By employing the RG technique, the 
present paper includes fluctuations past the mean-field BCS theory in a more sophisticated 
way than the original crossover theory developed by Nozieres and Schmitt-Rink S]. 
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critical reading of this paper. This work was supported by funds from Ministry of Education 
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